General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




mx3L~njm mm m.mam to smooc-ikjrtex iirrEiwcTioMs 

(HASA>TH-8541 1) SHOCK-FITTED EOLBfi H83-33850 

SOIOTIONS TO SHOCK-VOBTEX INTEBACTIOMS 
(NASA) 13 p HC A02/HF A01 CSC! 01A 

Due las 
G3/02 36471 


Nimiel E. Salas 
Tteoaas A. Zaiig 
M. Y<Msuff Hussalnl 


Report No. 82-24 

A prelfarinary version of this report by the same 
authors appeared as TM 84481, flay, 1982. 



August 16, 1982 

I 

WSTITUTE nm COliWER APPlianONS IN SCIENCE AND EN6INEERIN6 
NASA Langley Research Center, Hampton, Virginia 2366S 


Operated by the 



UR1VS61T1ES 


RESEARCH ASSOCIATION 


SHOCK-FITTED EULER SOLDTIONS TO SHOCK-VORTEX IHTERACTIONS 


Manuel D. Salas 
NASAt Langley Reeearck Center 

Thomas A* Zang 
College of William and Mary 

M* Tousuff Hussain! 

Inetitute for Computer Applioatione in Science and Engineering 


ABSTRACT 

The Interaction of a shock wave with a hot spot, a single vortex and a 
vortex street Is studied within the framework of the two-dimensional 
con^resslble Euler equations* The numerical results obtained by the 
pseudospectral method and the finite difference MacCormack method are 
compared. In both the methods the shock wave Is fitted as a boundary of the 
computational domain* 
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Introduction 


In their recent paper, Pao and Salas [1] presented a finite difference 
solution to the two-dimensional Euler equations governing the phenomena of 
shock wave Interaction «rlth an Isolated vortex* Their study emphasized the 
acoustic aspects of the problem* Zang, Hussalnl, and Bushnell [2] extended 
this numerical approach to the problem of turbulence amplification In shock 
wave Interactions* In this work It was necessary to resolve rather complex 
fine-scale structure In order to draw meaningful conclusions about the 
transient processes that dominate turbulence amplification* The present study 
Is a continuation of these efforts In two directions* First, to develop a 
highly accurate pseudo-spectral method capable of resolving the crucial small- 
scale structure on e relatively coarse grid* Second, to gain Insight Into the 
nonlinear dynamics of the transient processes Involved in the passage of a 
shock wave over a single vortex, a vortex street and a hot spot* 

Spectral methods have been demonstrated [3] , [4] to be powerful 
alternatives to finite difference methods for the numerical solution of smooth 
flows* Recently, the work of Gottlieb, Lustman, and Orszag [5] and of Zang 
and Hussalnl [6] have shown their applicability to simple one-dlmenslonal 
compressible flows with shocks* The present paper discusses a Chebyshev 
pseudo-spectral method that has produced reliable, accurate and efficient 
solutions to complex, two-dimensional flows with a strong shock* Although the 
main emphasis Is on the spectral technique, solutions to the governing 
equations obtained by the well known, second-order finite difference method 
originated by MacCormack are also given* This is the method that was used In 
[1] and [2]* The finite difference results for the present set of problems 
were calculated on a very fine grid and are used here for comparison with the 
solutions obtained with the spectral method* 



Statement of the Problem 


The physical problem that we model corresponds to an Infinite, Initially 
planar normal shock wave moving from left to right Into a downstream region 
containing a flow field representative of one or more vortices, or a hot 
spot* In order to model the Interaction of the shock wave with some given 
flow field ahead of It, It Is only necessary to compute the flow field 
upstream of the shock* The physical domain, therefore, need consist only of 
the region between some left boundary h(t). Judiciously chosen such that it 
will be far from the Interaction region, and the shock wave front Itself 
Xg(y,t)* It Is mapped onto the computational domain by the transformation, 

X . « - 

x^(y,t) - h(t) » 

V tanh(av) I OF POOR QUALITY 

^ " 2 » 

T - t* 


The computational domain Is thus (X,Y) e [0,1] x [0,1]. Note the 

stretching that has been used to handle the Infinite extent of the lateral 

coordinate y» If the relative shock Mach number M Is sufficiently high 

s 

(M^ > 2*08), the flow upstream of the shock remains supersonic* In this 
case, the left boundary corresponds to a supersonic Inflow, and all dependent 
variables can be prescribed on It* However, If the relative shock Mach number 
Is low, then radlatlon~type boundary conditions are used at the left 

boundary* On the right, the computational region Is bounded by the shock 

wave* Downstream of the shock the flow field Is given analytically* The flow 
field Immediately upstream of the shock, as well as the shape and velocity of 
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the shock, are evaluated such that the Ranklne-Hugonlot jump conditions and 
the compatibility condition reaching the shock wave from the upstream side are 
simultaneously satisfied* 

The unsteady, two-dimensional, compressible, Euler equations In the 
computational plane are written In the form, 

Qt + + BQy - 0 . 

where Q - [P,u,v,S]^ and 
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The natural logarithm of the pressure, the speed of sound, and the entropy are 
represented by P, a, and S, respectively, and y ts the ratio of specific 
heats* The velocity In the x and y directions are n and v, 
respectively* All variables are normalized with respect to reference 
conditions at downstream Infinity, as In [1]. The contravarlant velocity 
components are defined by 


U - + uX + vX 

t X y 

and 


V - + uY + vY * 

t X y 


Subscripts denote partial derivatives with respect to the Independent 
variables* 
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Let k denote the time level and let At be the time step Increment* 
The time discretization of eq* (12) is then as follows: 

Q - [1 - 

- V2 tQ^ + (I - AtDq] , 


where the spatial operator L represents an approximation to Ad^ + 

k 

the finite difference HacCormack method, the operators L and L are 
evaluated as two points forward and two points backward differences in the 
predictor (left) and corrector (right) levels, respectively. In the pseudo- 
spectral method studied here, the solution Q is first expanded as a double 
Chebyshev series. 


Q(X,T,I) - t t Q (T)t <t)T (II). 

p-0 q-0 P’ f ’ 

where 


5 - 2X - 1, 


T) - 2Y - I, 


and and are the Chebyshev polynomials of degrees p and q. 

derivatives appearing in the spatial operators are then evaluated as 


The 
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1 
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P’ s .4+1 *>’ 

aH^ odd 


and 



2 . 


c ■ 1, p > 0. 
P 



The Qy derivative la evaluated In a similar fashion* 

The evaluation of the shock wave shape and velocity followed the same 
procedure described In [1], except that In the spectral formulation, the 
derivatives that must be evaluated on the upstream side of the shock are 
expressed as Giebyshev expansions. At the left boundary, all variables were 
specified for supersonic Inflow. For the case of subsonic Inflow, the two 
velocity components and the entropy were specified, while the pressure was 
computed from a quasl-one-dlmenslonal characteristic as described In [8] . 

The pseudo-spectral method has a tendency to develop slowly growing 
oscillations. Because of the global nature of this method they are spread 
over the entire flow field rather than being confined to the vicinity of sharp 
gradients. The underlying smooth solution can be recovered by a variety of 
filtering techniques. The results presented here were obtained by applying a 
von Hann window filter (see [6J for details) every 160 time steps. Another 
practical consideration Is the explicit time-step restriction. The Chebyshev 
collocation points are clustered near the boundaries. Thus, smaller time- 
steps must be used In the pseudo-spectral calculations. 
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Perhaps the simplest Interaction to consider Is that of a planar shock 
wave with a ho. spot* as shown In Pig* 1* The flow field downstream of the 
shock wave situated at x ■ 0 at t " 0 Is taken as a quiescent field whose 
temperature distribution a. Is given by 

where for this case tc - 0.25, r • .125, Xq ■ 0.5, - 0.0, and - 3 at 

t * 0. The vortlclty distributions obtained by the finite dlffererce method 
and by the spectral method after the shock wave has passed over the hot spot 
are shown In Fig. 3. The finite difference solution presented here, and In 
all cases that follow, was obtained with 75 mesh points In the X direction and 
50 In T. The spectral solutions were all obtained with 33 collocation points 
In the X direction and 17 In Y. There Is very little difference between the 
two solutions. Both show the two counter-rotating vortices upstream of the 
shock, Which Is typical of this Interaction. See [2] for more details on the 
physics. 

Figure 2 shows the velocity field for a single vortex about to Interact 
with a shock wave of the same Initial strength as In the previous case. The 
downstream conditions here are obtained by assuming a constant density field, 
calculating the velocity from the stream function, 

<i» - log /r^ (x-xq)^ + * 

the pressure from Bernoulli's relation, and the temperature from the equation 
of state. For the case shown In Fig. 2, the circulation k « 2 and the 


softening scale r - 0.1. This model approaches an Idealized incompressible 
point vortex at large distances but Is much smoother near the center. Figure 
4 shows the resulting pressure field after the shcck wave has passed over the 
vortex. Overall, the results are qualitatively very similar, although the 
spectral method seems to be resolving the pressure field more accurately. See 
[1] for more details on the physics. 

Finally, Figs. 5 and 6 show the results for the Interaction with the 
Kerman vortex street that simulates the conditions of the experiment reported 
In [7]. For this case, the stream function representing the vortex is given 
by the difference of and (l'_ where 

log[cosh(™ + (y ± V 2 b)^) - cos(^ (x ±V 2 c))]. 


To match the experiment, the circulation, core radius, shock Mach number and 
vortex separation parameters were determined as ic ■■ 0.186, r > 0.1, - 
1.3, c ■ .33, and b ■ .048. For this calculation, the Inflow Mach number was 
subsonic and radiation boundary conditions were applied at the left 
boundary. The results shown In Fig. 6 are In agreement with the 
experimentally observed [7] longitudinal compression and lateral elongation of 
the vortex field after passage through the shock. The finite difference 
results are noticeably smoother than the spectral ones. However, It Is well 
known that the Idealized Karman vortex street Is unstable Zor all but one 
special raLlo of the horizontal and vertical separations. The downstream flow 
has the stable ratio. Thus, It Is likely that the upstream flow Is physically 
unstable. It may be that the spectral results have captured this phenomena, 
but more computations are needed to settle the Issue* 
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rig» I i Surface plot of entropy 
distribution for a hot apot about 
to interact vlth a Mach 3 shock 
wave* 


Fig. 2 ; Velocity field of a single vortex 
about to interact with a Mach 3 shock wave 
(solid curve)* The velocity vectors repre- 
sent perturbation from the mean values* 




Fig. 3 : Vorticlty fields at t • 0*20 computed by pseudo-spectral (left) and finite 

difference (right) methods for a hot spot after interaction with a Mach 3 shock 
wave (solid curves). 




Tit» 4 ; Isobars at t ■ 0.20 computed by pseudo-spectral (loft) and finite 
difference (right) methods for a single vorlcx after Inlriactlon with a 
V.^rh 3 shock wave. 
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Tt&« 5 t Velocity field (left) and vortlclty contoure (right) for a Karman vortex 
atreet about to interact with a Mach 1.3 ahock wave (solid curves). The velocity 
vectors represent perturbations from mean values. Negstlve contour levels sre 
are indicated by dashed lines. 



Flgj;__6; Veloclcy fields (top) and vortlclty contours (bottom) at f. - 0.36 computed 
by pseudo-spectral (left) and finite difference (right) methods for a Karman vortex 
street after interaction with ■ Mach 1.3 shock wave (solid curves). The velocity 
vectors represent perturbations from mean values. Negative contour levels sre indi- 
cated by dashed lines. 







